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We investigate relativistic and quantum electrodynamic effects for highly-excited bound states 
f— ^ ■ in hydrogenlike systems (Rydberg states). In particular, hydrogenic one-loop Bethe logarithms are 

calculated for all circular states (I = n — 1) in the range 20 < n < 60 and successfully compared 
£*0 _ to an existing asymptotic expansion for large principal quantum number n. We provide accurate 

expansions of the Bethe logarithm for large values of n, for S, P and circular Rydberg states. These 
three expansions are expected to give any Bethe logarithms for principal quantum number n > 20 to 
an accuracy of five to seven decimal digits, within the specified manifolds of atomic states. Within 
the numerical accuracy, the results constitute unified, general formulas for quantum electrodynamic 
corrections whose validity is not restricted to a single atomic state. The results are relevant for 
accurate predictions of radiative shifts of Rydberg states and for the description of the recently 
investigated laser-dressed Lamb shift, which is observable in a strong coherent- wave light field. 

PACS numbers: 12.20.Ds, 31.30. Jv, 31.15.-p, ll.10.Jj 

I. INTRODUCTION 

q , Circular Rydberg states (n — 1 = / = \m\) have attracted attention in the past two decades, in part because 

the transitions among these states are resonant with typical modes of microwave cavities. A convenient mechanism 
for the preparation of these states has been described in 0, 0, Q When coupled to the modes of a microwave 
cavity, the metastable Rydberg atoms undergo Rabi oscillations, leading to a direct verification of the concept of field 

CIh' quantization |J] . The coupling of the resonant cavity to the Rydberg atom can be used for the generation of entangled 

Einstein-Podolsky-Rosen pairs of atoms [a] . A further field of interest has been the theoretical simulation of atomic 

I ' non-dispersive wave packets, which can be realized with the "assistance" of a background laser field with a specific 

£> , "Kepler" frequency (see, e.g., [fJLllSl])- 

This partial list of applications is supplemented here by an investigation of relativistic and radiative effects in highly- 
excited hydrogenic states, without any consideration of the additional interaction with a cavity. The electron density 
^^ ■ of circular Rydberg states (n — 1 = I = |to|) resembles Kepler circular orbits (see Fig.^). However, the characteristic 
shape of wave function is restricted to the maximum-|m| subcomponent (see also Fig.^l. For states with nonmaximal 
\m\ < I, the shape of the electron density is manifestly different. For all circular Rydberg states, the wave function is 
spread out considerably in comparison to lower states, on a length scale of (r) ~ const, x n 2 asohr (see Ref. |jj), where 
r is the radial coordinate, aeohr is the Bohr radius and n is the (high) principal quantum number. Relativistic and 
radiative effects in Rydberg atoms allow for a simultaneous perspective on four different physical regimes: relativistic 
and QED radiative effects in the transition in the borderline region between classical and quantum physics. 

In addition to the aforementioned reasons for studying Rydberg atoms, the current pa per is also inspired by a 
recently completed analysis of radiative corrections for laser-dressed states [l(J, uM 03 • I n El LL2| > & n experimental 

Q_|. setup based on the hydrogen 1S-2Pj (j = l J2, 3 /2) transition is discussed in detail. A coherent-wave light source with 
the appropriate frequency has recently become available 13] . Because the 15 ground state is fully stable against 
• i— i , radiative (spontaneous) decay, the incoherent atomic fluorescence (in the case of the IS— 2Pj transition) is described 
/\ ■ to excellent accuracy by the Mollow spectrum |1J|. A further promising setup for studying the laser-dressed Lamb 
shift would be provided by a transition from the metastable 2S state to a high-n (Rydberg) P state. Rydberg 
states have a small radiative decay width, which would facilitate the observation of the predicted dynamic corrections 
to the Lamb shift [I(J,[l5|]. Furthermore, the frequency required for the excitation process 2S <^> nPj (high n) is 
smaller and therefore a priori less problematic to realize experimentally than the 1S-2Pj transition. In the formalism 
of [III UlL Ha], the self-energy correction analyzed here contributes to the Lamb shift of "bare" atomic levels; this 
contribution is an essential part of the laser-dressed Lamb shift, which may be observed in incoherent light scattering 
by hydrogenlike atoms. 

Another motivation for the current work is to find actual values of the Bethe logarithms for highly-excited states, 
since the Bethe lo gari thm gives a significant contribution to radiative level shifts. This is true of hydrogenlike, as 
well as heliumlike jig . J17L fl8| and lithiumlike [T]j systems. One of the most extensive systematic studies |2(j| of 
hydrogenic (one-loop) Bethe logarithms currently available in the literature extends only up to n = 20, and only 
approximate values are available for higher excited circular Rydberg states [23j. Furthermore, it is of interest to study 
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FIG. 1: (color online.) Plots of the radial probability density r 2 \ip(r, 9, <f>)\ 2 for two (nonrelativistic) states with quantum 
numbers n = 10 and / = 9, but different magnetic quantum numbers. Figure (a) displays the density of the states with 
\m\ — 9, in the plane of constant polar angle 9 = 7r/2. In space, the probability density of the electron in this maximum-|m| 
state describes a circular shape, which is why circular states are being viewed as analogues of classical planetary motion in 
the quantum domain. However, this characteristic circular pattern is restricted to the extremal magnetic quantum numbers 
(\m\ = I). Figure (b) shows the radial probability density of the state with m = 0, in the plane of constant azimuth = 0. 
The two maxima in the probability density are in the direction of the polar angles 9 — and 9 — n. 



the asymptotics of quantum electrodynamic corrections as a function of the quantum numbers. Such results can 
for instance be useful in comparing different calculational approaches to these corrections, and for obtaining general 
expressions valid for arbitrary quantum numbers within a specific manifold of states. Certain asymptotic structures 
have recently been conjectured for radiative corrections |2ll 1221 123 | . 

Furthermore, the analysis of highly-excited (bare) atomic states is motivated by recent experiments: e.g., in |24j . 
the evaluation of self energies of circular states of orbital quantum number I ~ 30 was required. The consideration 
of quantum electrodynamic corrections for highly-excited states is also relevant, in part, to the analysis of quantum 
electrodynamic corrections to dielectronic recombination resonances |'25j. 

In this work, we thus study radiative corrections for highly excited hydrogenic states. Specifically, after recalling 
basic facts about relativistic corrections in Sec. Ill Al we proceed to the study of radiative effects in Sees. Ill Bl and lllCl 
with an emphasis on the radiative decay width of circular states in Sec. Ill CI The main part of the investigation 
reported here is contained in III Dl where a numerical calculation of Bethe logarithms for highly and very highly 
excited circular states is described (see Table H] below). This investigation then leads to the asymptotic formulas for 
Bethe logarithms which yield a relative accuracy of 10 -5 or better, for general states within a specified manifold. 
The self-energy corrections studied in Sec. Ill Dl are the by far dominant radiative corrections for circular Rydberg 
states. Conclusions are drawn in Sec. lIIII This paper follows the usual convention for hydrogenic quantum numbers: 
the principal quantum number is denoted by n, the orbital angular momentum by I and the total electron angular 
momentum by j. As is customary in the literature, Z denotes the nuclear charge number of the hydrogenlike ion 
under consideration, and a is the fine-structure constant. Natural units (ft = c = eo = 1) are used throughout the 
text. 



II. RELATIVISTIC AND RADIATIVE EFFECTS 



A. Relativistic corrections 



We briefly recall some basic facts about relativistic corrections to highly-excited states in hydrogenlike systems. 
The numerical values of the coefficients multiplying the relativistic corrections of order (Za) 4 mc 2 are typically small 
for high-n states. In order to illustrate this point, we briefly recall the expansion of the Dirac bound-state energy in 
a hydrogenlike system in powers of Za up to order (Za) 4 (see, e.g., [2fJ, Eq. (2-87)]): 
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The (Za) 4 -term is the dominant relativistic correction for \ow-Z hydrogenlike systems and is responsible, in particular, 
for the fine structure. By inspection of Eq. JIJ, we infer that for circular states (which have j = n — 1 ± 1/2), the 



(Za) 4 -term scales as j~ 4 ~ n~ 4 , which contrasts with the familiar n _3 -scaling valid for fixed j. As we will see below, 
radiative corrections in Rydberg states also tend to display faster asymptotic behaviors (in n) than states with a fixed 
angular momentum. 

B. Self Energy and Radiative Lifetime 

In this article, we study radiative effects in highly-excited hydrogenic energy levels. Our focus is on the self-energy 
and the radiative lifetime. As pointed out in Ref. [27f, the self-energy and the radiative lifetime are intimately related 
to each other. Namely, the radiative lifetime of a hydrogenic state is proportional to the imaginary part of the 
ex pect ation value of that state {il)\Yi ren \^p) , where E r on is the renormalized self-energy operator discussed in Ch. 7 
of J23, and tp is the relativistic (Dirac) wave function. The self-energy is just the real part of (ip\Y, ren \ijj) and leads to 
the familiar shift of the energy levels. 

As is well known, the real part of the one-loop self-energy shift AE$e may be written as 

a i Z a I vrt 

AE SE = Re<V>|£ ren |V) = - l 3 F{nl h Za) . (2) 

Here, F(nlj, Za) is a dimensionless quantity and the the notation nlj follows the usual spectroscopic characterization 
of the hydrogenic state ip. It is customary in the literature to suppress the dependence of F on the quantum numbers 
n, j and I and write F{Za) for F(nlj, Za). 

The decay width T is the reciprocal of the lifetime r and is given by 

r=- = -2Im(^|S rcn K/;). (3) 

T 

For states with the same angular momentum quantum numbers, the decay rate decreases with n approximately as 
n~ 3 , and the lifetime correspondingly goes as n 3 (see Ref. [tj)- This is also manifest in the definition J5J of the scaled 
self-energy function F(nlj, Za), the structure of which reflects the usual, familiar scaling of the effect with n (for fixed 
I and j). 

C. Radiative Lifetime: Asymptotics for Rydberg States 

To a good approximation, the lifetime of a highly-excited hydrogenic state does not depend on the spin of the 
electron, and it is only the principal quantum number and the orbital angular momentum which enter into the 
leading-order (nonrelativistic) expression for the decay rate of excited hydrogenic states. 

We denote the (leading-order) decay width of a state with principal quantum number n and orbital angular mo- 
mentum I by the symbol T(n, I). Here, we are particularly interested in states with n = 1 + 1, and our goal is to study 
the dependence of the lifetime on the bound-state quantum numbers. For the manifold of states with n = I + 1, based 
on the formalism introduced in [27j . we obtain the following general expression, 

2a {Za) 4 m \ n^+\n~lf^ \ 
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The asymptotic expansion Ij4bfl about large n works surprisingly well, even at low principal quantum number. For 
n = 2, the first few terms listed in (|4b|) reproduce the exact formula (|4a|l to within 5 % accuracy. For n = 40, the 
difference is less than 10~ 10 in relative units. Equation (|4b|) implies that the lifetime of circular Rydberg states with 
I = n — 1 scales effectively as n 5 , in contrast to the familiar n 3 -scaling which is valid for fixed I, but varying n. 

Based on the well-behaved asymptotic structure of the leading non-relativistic contribution to lm('0|£ rcn |'!/>) for 
high quantum numbers, we now turn our attention to the real part Re(V'|S ron |-0). 



D. Self-Energy 

Vacuum-polarization effects are negligible for high-^ circular Rydberg states [the radial component of the wave 
function scales as (Zamr) 1 for small r\. The dominant radiative correction to hydrogenic energy levels is given by the 
(one-photon) self energy 
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which is a process in which the bound electron (double line) emits and re-absorbs a photon (wavy line); this effect 
shifts the energies predicted by the Dirac equation. 

We now turn to the investigation of the energy shift A-Ese, which is defined in Eq. J5J and corresponds to the 
above diagram. The semi-analytic expansion of F(nlj, Za) about Za = for a general atomic state with quantum 
numbers n, I and j gives rise to the expression 28] 



F(nlj, Za) = A 41 (nlj) hx[(Za)~ 2 ] + A iQ {nl ) + O(Za) . 



(5) 



The Za-expansion is semi-analytic, i.e., it involves powers of Za and of ln[(Za)~ 2 ]. The A coefficients have two 
indices, the first of which denotes the power of Za [including those powers contained in Eq. J2Jl], while the second 
index denotes the power of the logarithm ln[(Za)~ 2 ]. 
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FIG. 2: (color online.) Relative difference between the Bethe logarithm values presented in Tab. U] and the asymptotic 
formula ©, for circular Rydberg states (I — n — 1) in the range 20 < n < 60. Since zero is contained in all the error bars, 
the Bethe logarithms in Tab. [I] are fully consistent with the approximation ©. The error bars in the figure correspond to the 
uncertainty in the numerical coefficients of Eq. Q. 

The coefficient A^ jnlj] assumes a value of 4/3 for S states and vanishes for all non-S" states. The general formula 
for A i0 reads (see, e.g., [29tl30| l 
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where k = 2 (l—j) (j + 1/2). The Bethe logarithm hiko(n, I) is an inherently nonrelativistic quantity, whose expression 
in natural units reads 
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Here, Hg is the nonrelativistic Coulomb Schrodinger Hamiltonian, p l is the ith Cartesian component of the momentum 
operator, E n and <f> are respectively the (nonrelativistic) energy and the wave-function of a state with the quantum 
numbers [n, I), and the summation over Cartesian coordinates [i — 1, 2, 3) is implicit. 

Bcthe logarithms have been studied for all hydrogenic states with n < 20, by various authors and with increasing 
accuracy due to advances in the algorithms used and in the computing technology available [2(J, l3ll l32l l33l J34L l35l 13a . 
|^3j|2a|. While the calculational difhculties (numerical convergence properties of the involved hypergeometric functions) 
increase as the principal quantum number increases, the evaluation of the hydrogenic nonrelativistic Bethe logarithm 
can be regarded as an easy computational task from the perspective of our current understanding of bound states, and 
of modern computer technology. In the current work, we employ the technique previously discussed in [23|, |3!j , where 
the relevant matrix element of the hydrogenic wavefunction |4l| is expressed in terms of hypergeometric functions 2 -Pi) 
which arc then suitably integrated. For high energies of the virtual photon, we found it useful to employ the combined 
nonlinear-condensation transformation |40J, |42j in order to accelerate the convergence of the series representation of 
these hypergeometric functions. 

One of the results we have obtained reads 

lnk (n = 20,l= 19) = -0.000 008 084 977 837 087 891(1), (8) 

a value which confirms the result obtained previously in Ref. |2(J] for this state. The error is due to the uncertainty in 
the final integration over the virtual photon energy. While the numerical accuracy of this result is merely of academic 
interest, we found it useful to verify the accuracy of our computational method against the basis-set techniques 
employed in |20ll3q . 



TABLE I: Values of the Bethe logarithm for circular Rydberg states (for which 
I = n — 1), in the range 20 < n < 60. According to Eqs. © and ©, the Bethe 
logarithm lnfeo contributes significantly to the radiative correction to the energy of 
hydrogenic states. All explicit results found here agree with a previously published 
expansion |23l Eq. (37)], as well as with the new asymptotic expansion ©, as is 
shown in Fig. (within the numerical uncertainties in the asymptotic coefficients). 
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-0.69410 65660 642(1) 
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-0.08782 81749 454(1) 


x 10 _b 


22 


-0.60031 18198 591(1) 


x 10~ 5 


42 


-0.08158 06567 501(1) 


x 10~ 5 


23 


-0.52267 98437 374(1) 


x 10~ 5 


43 


-0.07591 18687 638(1) 


x 10~ 5 


24 


-0.45787 44770 701(1) 


x 10~ 5 


44 


-0.07075 63248 884(1) 


x 10" 5 


25 


-0.40335 28104 402(1) 


x 10~ 5 


45 


-0.06605 72331 309(1) 


x 10~ 5 


26 


-0.35715 36036 972(1) 


x 10 -5 


46 


-0.06176 51782 987(1) 


x 10~ 5 


27 


-0.31774 78119 539(1) 


x 10~ 5 


47 


-0.05783 70279 429(1) 


x 10 -5 


28 


-0.28393 14416 323(1) 


x 10~ 5 


48 


-0.05423 50198 720(1) 


x 10~ 5 


29 


-0.25474 77154 839(1) 


x 10" 5 


49 


-0.05092 59980 293(1) 


x 10~ 5 


30 


-0.22942 98276 714(1) 


x 10~ 5 


50 


-0.04788 07701 395(1) 


x 10~ 5 


31 


-0.20735 83589 839(1) 


x 10~ 5 


51 


-0.04507 35657 214(1) 


x 10 -5 


32 


-0.18802 92605 357(1) 


x 10 -5 


52 


-0.04248 15771 656(1) 


x 10~ 5 


33 


-0.17102 95447 205(1) 


x 10" 5 


53 


-0.04008 45698 251(1) 


x 10" 5 


34 


-0.15601 86575 648(1) 


x 10~ 5 


54 


-0.03786 45496 597(1) 


x 10~ 5 


35 


-0.14271 40813 819(1) 


x 10~ 5 


55 


-0.03580 54790 517(1) 


x 10 -5 


36 


-0.13088 01169 928(1) 


x 10 -5 


56 


-0.03389 30330 782(1) 


x 10~ 5 


37 


-0.12031 90769 194(1) 


x 10 -5 


57 


-0.03211 43898 736(1) 


x 10~ 5 


38 


-0.11086 43219 585(1) 


x 10 -5 


58 


-0.03045 80498 130(1) 


x 10~ 5 


39 


-0.10237 47182 006(1) 


x 10~ 5 


59 


-0.02891 36791 317(1) 


x 10" 5 


40 


-0.09473 01966 836(1) 


x HT 5 


60 


-0.02747 19743 302(1) 


x 10" 5 



ar Rydberg states, for 21 < 
in Ref. [H Eq. (37)]. This 
where the range 1 < n < 20 



In Tab. H] we present accurate numerical values for the Bethe logarithm Q of circu 
11 < 60. These values are in agreement with the truncated asymptotic expansion found 
asymptotic expansion had been derived on the basis of numerical data reported in Ref. [20( 

had been covered. Based on the numerical data of Tab. [IJ which cover a wider range of principal quantum numbers 
as compared to Ref. [2(| , we obtain the following improved asymptotics for the Bethe logarithm of circular Rydberg 
states: 

0.0248208(6) 0.03814(2) 0.1145(5) 0.166(3) 0.22(2) 
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where terms of order l~ k with k > 6 are neglected. The algorithm described in the Appendix of [23 was used 
for obtaining the coefficients of this asymptotic expansion. All the coefficients in © separately are in very good 
agreement with those found previously J23j Eq. (37)]. The Bethe logarithm values in Tab. Q] are fully consistent with 
the above truncated expansion: they fall inside the error bars of ©, as is illustrated in Fig. [2 Furthermore, the 
"expectation" values for the coefficients of formula © correctly "predict" about seven digits of the actual numerical 
value of the Bethe logarithm in Tab.HJ as illustrated in Fig.|3 Based on the very good consistency of the asymptotic 
expansion with the actual numerical data (Fig. |2J), we conjecture that the "expectation" values of the expansion @ 
give the Bethe logarithm of all circular Rydberg states with n > 20 with a relative precision better than 10 -6 , as can 
be expected from Fig. [5J Furthermore, we would like to conjecture here that the expansion (|5J not only represents a 
"polynomial fit" to the numerical data, but that it represents the true asymptotic expansion of the Bethe logarithm 
for high quantum numbers. 

For S states, based on data available for principal quantum numbers n < 20 (see Ref. |2(J), we obtain the following 
asymptotic expansion as a function of the principal quantum number n, 

ln fc „(V = 0) * 2.72265434(5) + °- 000000(5) + ^f^ ~ ^"^ + ^ ~ ^ , (10) 

n n z n 6 n 4 n b 

where terms of order n~ k with k > 6 are neglected. By a reasoning similar to the one above, we expect the 
expansion JTDJ, with the "expectation" values employed for the coefficients, to give the Bethe logarithm of all S states 
with n > 10 to better than 2 • 10~ 7 in relative units. 
Finally, for high-n P states, we obtain 

iw 1 n nruon^m , P- 000000 ^) , 0-20530(15) 0.599(5) 1.45(10) 3(1) 

lnfc (n, 1 = 1) ~ -0.0490545(1 H 1 5 5 1 -, =- . (11) 

n n z n 6 n 4 n° 

where again terms of order n~ k with k > 6 are neglected. We expect Eq. Ill 111 to reproduce the exact Bethe Logarithms 
of all P states for n > 10 to a relative accuracy better than 5 • 10 -5 . 

III. CONCLUSIONS 

We have studied relativistic and radiative corrections for highly-excited hydrogenic states. In addition to S and 
P states, we have concentrated on circular Rydberg states (for which I = n— 1). We have observed on the relativistic 
and radiative corrections considered here that circular Rydberg state corrections have a faster asymptotic behavior 
(as a function of n) than states with a fixed angular momentum. 

The following scaling properties hold for the dominant relativistic and radiative effects on circular hydrogenic 
Rydberg states: (i) The leading relativistic corrections to the energy levels (of order (Za) 4 mc 2 ) scale as n~ 4 . (ii) 
Radiative decay rates scale as n~ 5 [see Eq. (I4b|l ]. (Hi) The Bethe logarithm contribution to the self-energy shift (J2J) 
scales as n~ 6 . Note the three inverse powers of I in Eq. @ which have to be taken into consideration in addition to 
the three inverse powers of n originating from @. By contrast, the Bethe logarithm contribution for highly-excited 
S and P states scales as n~ 3 [see Eqs. (|10|) and Ijllj l]. (iv) The anomalous magnetic moment contribution to the 
Lamb shift [which corresponds to the term — 2 , *, , in Eq. ©] scales as n~ 5 for circular Rydberg states and thus 
dominates over the Bethe logarithm in the limit of large I = n — 1. 

In this article, we also provide an exact result for the decay rate of highly-excited hydrogenic Rydberg states [see 
Eq. (|4a|) ]. Furthermore, we obtain accurate values (see Tab.QJ for the Bethe logarithms of circular Rydberg states (for 
which the orbital angular quantum number I is related to the principal quantum number as I = n— 1). These values 
are in a range of quantum numbers 21 < n < 60, where no such data have been available up to now, to the best of our 
knowledge. These results confirm that the dependence of radiative corrections on the bound-state quantum numbers 
can be represented, to a very good accuracy, by asymptotic expansions about large n which involve negative integer 
powers [sec Eqs. Il4b|l . (EJ, Ijl0|l and (|1 1 f> ] - This functional dependence has recently been observed for a number of 
QED effects |2lJ, |22j |23j . We conjecture here that the Bethe logarithm expansions ©, i|10|) and (|11|) represent the 
true asymptotic expansions about large quantum numbers (to the precision indicated in the text). 

The results obtained here are relevant for any analysis (e.g., J24|) where the evaluation of the self energy of circular 
(n = I + 1) states of orbital quantum number Z > 20 is required. They are also relevant for the description of the 
"bare" Lamb shift contribution to the Lamb shift of laser-dressed states [T(j, HH [H, Ej] in experiments with Rydberg 
states. A possible experimental setup would involve the metastable 2S state and an excitation to a high-n P state. 
The required frequency is smaller than the one required for the excitation IS <^> 2Pj, which has been studied in 
detail in [IJ. Moreover, both states involved in the transition 2S 43> nPj (high n) have a long lifetime, which would 
facilitate the observation of the dressed-state corrections to the Lamb shift [l(J, fill Il2| . The formalism required for 



the analysis of further corrections to the laser-dressed Lamb shift has recently been described in |l2( . The additional 
corrections are due to ionization, counter-rotating terms (Bloch-Siegert shifts, see [43J1 and nonresonant levels, as 
well as laser-field configuration dependent effects. 

Finally, we note that an accurate investigation of the Lamb shift of the (numerous) different magnetic subcom- 
ponents of circular Rydberg levels would allow for a sensitive test of the recently proposed, somewhat speculative 
"polarized Lamb shift" (see Ref. J44J) that predicts a breaking of the rotational symmetry imposed on quantum 
electrodynamics. 
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